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ABSTRACT 

New generation TeV gamma-ray telescopes have discovered many new sources, 
including several enigmatic unidentified TeV objects. HESS J0632+057 is a par- 
ticularly interesting unidentified TeV source since: it is a point source, it has a 
possible hard-spectrum X-ray counterpart and a positionally consistent Be star, 
it has evidence of long-term VHE flux variability, and it is postulated to be a 
newly detected TeV/X-ray binary. We have obtained Swift X-ray telescope ob- 
servations of this source from MJD 54857 to 54965, in an attempt to ascertain 
its nature and to investigate the hypothesis that it's a previously unknown X- 
ray/TeV binary. Variability and spectral properties similar to those of the other 
3 known X-ray/TeV binaries have been observed, with measured flux increases by 
factors of ~3. X-ray variability is present on multiple timescales including days 
to months; however, no clear signature of periodicity is present on the timescales 
probed by these data. If binary modulation is present and dominating the mea- 
sured variability, then the period of the orbit is likely to be > 54 days (half of this 
campaign), or it has a shorter period with a variable degree of flux modulation 
on successive high states. If the two high states measured to date are due to 
binary modulation, then the favored period is approximately 35-40 days. More 
observations are required to determine if this object is truly a binary system and 
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to determine the extent that the measured variability is due to inter-orbit flaring 
effects or periodic binary modulation. 

Subject headings: gamma rays: observations — X-rays: binaries — X-rays: indi- 
vidual(HESS J0632+057) 



Introduction 



During the past 5 years, TeV gamma ray astronomy has experienced rapid growth due 
to the commissioning of a new generation of sensitive atmospheric Cherenkov telescopes. 
These telescopes have led to a n expansion of the TeV catalog from ~10 sources to nearly 



100 sources in just a fe w years (IHoran fc Wakely 1120091 ; IWeekes 1120091 ; iBuckley et al. 112008 



Aharonian et al. II2008I ). Many of these sources have come from the survey of the Galactic 
plane by the HESS telescopes in Namibia. One of the most exciting results of this survey is 
a large number of unidentified TeV sources, which offer a new region of discovery space at 
very high energies. HESS J0632+057, located within the Monoceros Loop region, is one of 
these unidentified sources. 

HESS J0632+057 was first detected by HESS using 13.5 hours of data collected between 
March 2004 and March 2006 (lAharonian et al. 1 120071 ). It has a peak significance in the 
field of view of 7.1a, and the post-trials significance is 5.3a, after accounting for the trials 
associated with searching the entire field of view. The spatial distribution is consistent 
with a point source (RMS size < 2 arcmin at 95% confidence) near the edge of the Rosette 
Nebula, located at RA 06h 32' 58.3", Dec. +05° 48' 20" (± 28" stat., 20" sys.). The source 
position is located on the edge of the 99% contour of the EGRET unidentified object 3EG 
J0634+0521, but no object at the p osition of HESS J0632+057 is listed in the the Fermi 



bright so urce list (lAbdo et a 



2009). MWC 148, a Be star, is at the centroid of the HESS 



position. I Aharonian et al. I (120071 ) suggested that the TeV gamma ray emission could be 



associated with either: a) a shock front driven by the stellar wind of MWC 148, b) shock 
accelerated cosmic rays interacting with a molecular cloud at the edge of the Monoceros 
Loop, c) an association with an unidentified ROSAT source 1RXS J063258. 3+054857, d) 
an association with EGRET unidentified gamma ray source 3EG J0634+0521, and/or e) a 
previously unknown X-ray binary source associated with MWC 148. 

On 2007 September 17, X MM- Newton observed the region surrounding HESS J0632+057, 
obtaining 26 ksec of good data (IHinton et al. 1120091 ) . This resulted in the detection of several 
sources, the brightest being a point source, XMMU J063259. 3+054801, positionally coinci- 
dent with MWC 148 and HESS J0632+057 and was within the 99% error circle of 1RXS 
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J063258. 3+054857. During the relatively short observation period, the object exhibited sig- 
nificant variability in the form of a gradual decline from ~130 c/ksec to ~90 c/ksec. During 
this observation, it had a mean deabsorbed 1-10 keV flux of (5.3 ± 0.4) xl0~ 13 erg cm -2 
s _1 . The spectrum was well fit with an absorbed power law spectrum, and it had a hard 
photon index of 1.26 ± 0.04, similar to the spectra of known TeV binary systems. 

Follow-up observations of HESS J0632+057 at gamma-ray energies above 1 TeV by 
VERITAS during December 2006 - Jan 2007 and Dece mber 2008 - January 2009 yielded 



upper limits well below the fluxes published by HESS ( lAcciari et al. Il2009bl ). The non- 
detection of HESS J0632+057 by VERITAS together with the HESS results provides evidence 
for variability in the gamma-ray flux on time scales of months, although the available data 
do not allow any conclusion on a possible periodicity. 

At optical wavelengths, no reports of variability have been found in the literature. 
Significant flux variabili ty is reported in rad io at 5 GHz on about month time-scales around 
a mean flux of 0.3 mJy (ISkilton et al. Il2009i ). but no periodic variability could be discerned 
from these data. This radio flux is orders of magnitude lower than the typical radio flux 
expected from a TeV blazar which makes that potential interpretation unlikely. 

The point-like nature of the detected TeV source, the e xcellent positiona l coinc idence 



with MWC 148 (chance coincidence of ~ 10 4 according to lAharonian et al. I (120071 )). the 



location on the Galactic plane with a low radio flux, and the variable X-ray and gamma- 
ray emission from a location coincident with MWC 148 are all facts that argue in favor 
of an X-ray binary in association with MWC 148. There are currently just three known 
binary systems that are confirmed as pro ducers of de tected TeV gamma ray e mission, namely 
PSR B1259-63 (kharonian et al. Il2005bh. LS 5039 JAharonian et al. Il2005ah . and LS I +61 



303 (Albert et al. 



2006 



Acciari et al. 



20081 . |2009i ). There is also an accreting black hole 



binary system Cyg nus X-l, with evidence of detection of a GeV/TeV gamma-ray flare 
( {Albert et al. 1120071 ). but this result relies upon marginal significance which would benefit 
from future confirmation if a high state occurs in the future. The three confirmed TeV 
binaries are powered by pulsar winds driving shock acceleration or by accretion onto the 
compact object driving a microquasar jet. Confirmation of HESS J0632+057/MWC 148 
as a TeV binary would add a 4th object to this short and special list and allow us to 
begin studying its interesting properties; whereas, refutation of the binary hypothesis would 
establish this as an even more mysterious and interesting class of TeV source. The most 
direct way to test this binary hypothesis is to search for periodic emission signatures. In this 
paper, we report on recent monitoring data taken with Swift-XRT with the aim of advancing 
our understanding of this source. 
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The Observations 



Swift was launched on 2004 November 20 (IGehrels et al. 1 120041 ) with a primary goal 
to provide detailed measurements of numerous gamma ray bursts (GRBs) and their after- 
glows with unprecedented reaction times. An ability to slew quickly and frequently has also 
made Swift an excellent multiwavelength observatory for many non-GRB sources through 
pre-planned monitoring observations and rapid response target of opportunity (ToO) ob- 
servati ons. The narrow fiel d instruments used for these studies are the X-ray telescope 



(XRT ; iBurrows et al. I (120051 )) and the Ultraviolet-Optical Telescope (UVOT; iRoming et al. 
( 120051 )). Bright, hard-spectrum sources can also be studied with the burst alert telescope 
(BAT) in the ~ 15 - 150 keV energy band. For these HESS J0632+057 observations, UVOT 
was overwhelmed by the bright star in the field of view. The source was not bright enough 
to be detected by BAT at high energies, or to allow BAT to produce a constraining upper 
limit. All of the analysis in this paper is based upon XRT data. 

The XRT observations were typically taken as 4-6 ksec exposures. The initial set of ex- 
posures, from MJD 54857 to MJD 54874, were taken with various sampling intervals chosen 
primarily to overlap with th e simultaneous observ ations at the VERITAS very high energy 



gamma ray observatory (see lAcciari et al. I (j2009bl ) for results) . These observations included 



a series of 4 observations spaced on daily time intervals, as well as two more observations 
with 1 week spacing. One month later, starting on MJD 54904, a series of monitoring obser- 
vations was initiated with ~3-4 day time spacing. This resulted in 20 additional observations 
extending from MJD 54904 to MJD 54965, at which time the spacecraft had to discontinue 
observations due to a Sun constraint. The total Swift-XKT data set includes 103 ksec of 
observations extending from 2009 January 26 to 2009 May 20. All observations beyond 
January 29 are being reported here for the first time. 



3. Analysis 

The Swift-XBT data were processed using the most recent versions of the standard 
Swift tools and the most recent calibration files available at the time of data processing. 
In particular, we utilized Swift Software version 3.0, FTOOLS version 6.6.3, and XSPEC 
version 12.5.0. Light curves were generated using xrtgrblc version 1.3. 

All of the observations were obtained in photon counting (PC) mode. Circular and 
annular regions are used to describe the source and background areas respectively, and the 
radii of both regions depend on the measured count rate. In order to handle cases where 
the sources land on bad CCD detector columns, point spread function correction is handled 
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using xrtlccorr. Since these observations always resulted in XRT count rates ~0. 01-0. 06 c/s, 
there was no significant pile-up, which occurs in photon counting mode at count rates above 
~0.5 c/s. The full light curves use a bin size of one observation per bin (observations were 
typically 4-5 ksec duration). We also evaluated single observation light curves, and in these 
cases, we used a bin size minimum of 20 cts/bin. All error bars are reported at the 1-sigma 
level, unless otherwise specified. 

As described in the results section below, we calculated spectral fits for the entire com- 
bined data set. We also calculated spectral fits for a series of observations that were binned 
temporally with large enough bins to provide enough photon statistics to allow reasonable 
spectral fits. Following the temporal binning, each observation was fit spectrally using a 
minimum energy binning ratio of 12 photons/bin, and in most cases the binning ratio was 
20 photons/bin. Once the spectral fits were obtained, the flux was calculated. These flux 
values, which were calculated for each of the temporally binned data segments, were then 
used to calculate a flux-to-rate ratio during that time period. This ratio was then applied to 
the rate light curve to convert it to a flux light curve, as shown in Figure [IJ which assumes 
that any spectral variability has a negligible impact on the flux-to-rate ratio within a single 
time bin. 



4. Results 

The X-ray light curve is shown in Figure [TJ There is clear evidence of strong variability 
on multiple timescales, with measurable flux doubling on timescales as short as ~5 days. 
While shorter timescale variability is certainly possible (within the individual ~5 ksec ob- 
servations), it is not measurable within the flux limited error bars provided by these data. 
Longer timescale variability is also evident, with an extended rise and decay lasting for at 
least 30 days within the time frame bracketed by MJD 54907 and MJD 54941. 

For this work, the hardness ratio is defined as R2/R1, where R2 is the rate in the 2-10 
keV band and R\ is the rate in the 0.3-2 keV band. The hardness ratio as a function of 
time is shown in Figure [2j The mean hardness ratio is 0.96 with a standard deviation of 
0.27. Fitting the hardness curve to a constant in time results in a reduced x 2 equal to 2.01 
with 24 degrees of freedom, with a chance probability of 99.8%. While the constant hardness 
hypothesis does not provide a good fit, it is not poor enough to claim significant evidence of 
spectral variability. 

A simple absorbed power law model, with Nh left to vary freely, was fit to the time 
binned data. The Nh was not found to vary significantly. The mean Nh was found to be 
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3.3 ±0.2 x 10 21 cm 2 , and all values are consistent, within error, with the value obtained from 
earlier XMM observations, 3.1 ± 0.3 x 10 21 cm -2 (IHinton et al. II2009I ). For all subsequent 



fitting, the Nh was fixed to the value obtained from the XMM observations. The photon 
indices from these fits are shown in Figure [3J The time axis error bars represent the full 
extent of the time bin used for each model fit, and the vertical error bars are 90% confidence 
level error bars. We also binned all of the data together to obtain a combined fit of an 
absorbed power law to the entire data set. The photon index was found to be 1.56±0.06, 
and the fit of the absorbed power law resulted in a reduced x 2 °f 0.89 with 72 degrees of 
freedom. The model and data for the overall fit are shown in Figure HI 

A power spectrum periodicity search was performed on the overall light curve shown 
in Fi gure [Tj using a Lo mb normalized periodogram which is applicable to unevenly spaced 
data (IPress et al. 1120071 ). However, no significant periodic signal can be discerned from these 
data. 



Discussion & Conclusions 



The X-ray flux from the TeV binary LS I +61 303 (which is potentially a very similar 
object to HESS J0632+057) varies over a 26.5 day orbital cycle at many w avelengths, in 



eluding variations by factors of 2 to 6 at both optical and X- ray wavelengths (IParedes et al. 



997 



Hol der. Falcone, fc Morris 



(ISmith et al. 



2009 



2007 ; ICasares et al. II2005I ). and it has fast timescale flares 



Barthelmy et al. 1 120081 ) . Therefore, we have probed HESS J0632+057 



for binary orbit variability timescales ranging from days to weeks, as well as searching for 
any large amplitude flaring on short timescales during all observations. The detection of a 
binary period would clearly allo w us to differentiate b etween X-ray binary models and Be 
star wind-driven shock models (ITownsend et al. 1 120071 ). as well as other less likely models 
such as cosmic rays interacting with a molecular cloud or a remarkably low radio flux blazar 
on the Galactic plane. 

While it is clear that further monitoring observations are needed in order to obtain any 
conclusions regarding periodicity, it is also clear that the X-ray light curve of this point 
source is variable on long (> weeks) timescales. The light curve is consistent with periodic 
variability on timescales longer than those probed in this study, if we assume that there 
is some additional variability due to flaring, but longer baseline observations are needed to 
evaluate the reality of any periodic variability. All of the measured fluxes in the epochs during 
which Swift observed the source are larger (by a factor of ~two) than the flux measured by 
XMM in September of 2007. A constant flux hypothesis is ruled out. A simplistic sinusoidal 
curve is also ruled out. However, it is possible that a more complicated model, involving 
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Fig. 1.— The X-ray light curve of XMMU J063259.3+054801 from the Swift-XRT observa- 
tions in the 0.3-10 keV band. The bottom panel shows the flux derived from binning up all 
photon events into time bins that are illustrated by the horizontal error bars and finding the 
best fit spectral model with N# fixed to 3.1 xlO 21 . The flux in the top panel explores tem- 
poral variability in more detail by deriving the flux directly from the X-ray rate light curve, 
with a flux to rate ratio derived for each of the spectral time bins defined by the horizontal 
error bars in the bottom panel (this assumes that spectral variability is a negligible effect 
within each time bin). 
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Fig. 2. — The X-ray hardness curve of XMMU J063259. 3+054801, where hardness is defined 
as the flux in the 2-10 keV band divided by flux in the 0.3-2 keV band. The mean hardness 
is shown as a solid line. A fit to this line produces a reduced y 2 of 2.01 with 24 degrees of 
freedom. 
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Fig. 3. — Spectral parameters from application of an absorbed power law model to time 
binned Swift-XKT observations of XMMU J063259.3+054801. The top panel is the power 
law photon index. The bottom panel is the reduced x 2 °f the model fit. Vertical error bars 
represent 90% confidence range. Horizontal error bars represent the extent of the time region 
used in the binning, which is the same temporal binning that was shown in Figure 1. 
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Fig. 4. — Absorbed power law model fit to sum of all Swift-XRT data analyzed in this paper. 
N H is fixed to the value obtained previously by XMM, and the photon index is found to be 
1.56±0.06, with a reduced x 2 of 0.89. Residuals are shown in the bottom panel. 



- 11 - 



a long period sinusoidal curve plus flares, could be fit to the light curve. The results of a 
Lomb-Scargle periodogram analysis show no significant periodicity; however, this analysis 
does not rule out the possibility of a periodicity at longer timescales beyond the range 
of these observations (>54 days) or the possibility of a shorter periodicity timescale with 
variable or small fractional flux modulation. For example, if the modulation fraction was 
<30%, then Swift-XKT would not be sensitive enough to see a significant deviation from one 
5 ksec pointing to the next, given the mean flux from HESS J0632+057. It is also likely that 
this study would miss periodic flux modulation fractions somewhat larger than this, due to 
the contamination of the periodic signal by non-periodic flare-like emission. Although the 
following is not statistically significant, we should also mention that if the two high states 
measured in these data are due to a binary modulation, then the favored period would be 
approximately 35-40 days; but, of course, there is no evidence for or against the periodic 
nature of these two high states since longer timescale observations are necessary for making 
such a determination. If these two high states are not related to any potential periodicity, 
then the favored period would be longer than half the timeframe of observations performed 
thus far (i.e. longer than 54 days). 

By plotting Nh/Nh^i versus binary orbit for known high mass X-ray binaries, Be- 
HMXBs and Sg-HMXBs have been found to occupy different regi ons of paramete r spac e 



(IBodaghee et al. 1120071 ; iRodriguez fc Bodaghee 1120081 ). According to lKalberla et al. I (120051 ) 



the galactic neutral hydrogen column density, Nh^u is 5.6 xlO 21 cm 2 in the direction of 
HESS J0632+057. If we take the orbital period of HESS J0632+057 to be >54 days (half 



the t ime period of these observations) and the Nh to be 3.1 x 10 cm (IHinton et al. 



20091 ). then it falls within the middle of the Be-HMXB section of the Bodaghee Diagram. 
This period, as well as the potential ~35-40 day period would both be long enough to be 
consistent with the interpretation of this object as a TeV emitting HMXB associated with 
the Be star MWC 148. 

If MWC 148 was a single star with a wind driven shock, it would be expected to 
have a reasonable spectral fit with a thermal model such as an absorbe d MeKaL model or 



an absorbed Raymond -Smith model with a temperature of kT~l keV (IRauw et al. 1 12002 



Townsend et al. 1120071 ). However, we find that such models do not fit the data well due to the 
presence of a harder component. The absorbed Mekal model results in a high temperature of 
kT=17±4 keV, which is too high to be likely for a single star model, with reduced x 2= l-36 for 
109 degrees of freedom. The absorbed Raymond-Smith model results in a reduced x 2= 8.15 
for 109 degrees of freedom. The spectrum is more naturally fit by an absorbed power law 
spectrum, which results in a reduced \ 2 °f 0-89 an d a power law photon index of 1.56±0.06. 
While not definitive, this evidence favors models involving more than a single star, such as 
the binary hypothesis. 
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These monitoring observations also provide the possibility to search for short timescale 
flaring from either XMMU J063259. 3+054801 or from any of the other dimmer X-ray sources 
in the field of view. X-ray flaring has been observed from other X-ray binaries and from at 
least one TeV emitting X-ray binary, namely LS I +61 303 (ISmith et al. 1 120091 ) . Short 
timescale flaring would provide a way to probe the size of the emission region, as well as the 
power of the engine that must be feeding the associated acceleration site. In principle, short 
timescale variability could also be used to probe different temporal signatures expected from 
Be star wind-driven shock models, binary interaction region models, and accretion- driven 
microquasar jet models. Variability as short as ~ 5 days was observed, but no flaring at 
shorter timescales was observed in the flux range probed by these observations. 

These data are complemented by the sim ultaneous data that h ave been obtained by 
Fermi, as we ll as radio obse r vations reported by ISkilton et al. I (120091 ) and TeV observations 
reported by lAcciari et al. I (J2009bJ). In future work, modeling of this composite spectral 
energy distribution and longer timescale observations will provide strong constraints on po- 
tential periodicities and on the acceleration site parameters, such as the B field and the 
electron distribution, implied by fitting a synchrotron spectrum plus an inverse Compton 
spectrum to the observations. 



This work is supported at Pennsylvania State University by NASA contract NAS5-00136 
and grant NNX08AV77G. 
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